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(57) Abstract: A method and an apparatus for extracting signals which are indicative of the level of consciousness of a patient 
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£>• model with exogenous input, and then calculating an index (AAI), which is displayed or used otherwise, indicative of the anaesthetic 
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Monitoring Auditory Evoked Potentials 
Background of the invention 

The present invention relates to a method and an apparatus for calculating 
an index representing the depth of anaesthesia when monitoring the audio 
evoked potentials of a person, and a computer program adapted to perform 
the method on a computer. 

Assessment of depth of anaesthesia is in general based on clinical obser- 
vations of physiologic parameters such as blood pressure, heart beat rate, 
pupil size etc. The use of neuro-muscular blocking agents during general 
anaesthesia disables the clinical signs that normally indicate conscious- 
ness. A number of incidents exist where patients describe that they were 
fully conscious during the surgery, and in the worst case had perception of 
pain and cardiac arrhythmias. Hence, there is a need for a method and ap- 
paratus to assess the anaesthetic depth. A number of investigation results 
have already been published, where Auditory Evoked Potentials (AEP) are 
used to indicate the level of consciousness during general anaesthesia. 
The AEP is a sub-component of the EEG signal, and it is elicited by acous- 
tic stimuli and is recorded with scalp electrodes, amplified and analysed by 
a computer. The AEP is an electrical, small signal embedded in noise from 
the ongoing EEG, and for this reason advanced signal processing is neces- 
sary to extract the AEP signals. The AEP signals are traditionally extracted 
by the averaging of up to 1000 repetitions of the response of the stimuli 
signals. This is a very time-consuming process, which takes up to several 
minutes to carry out, typically 2-3 minutes, which is excessive if the anaes- 
thesiologist has to use the AEP signals as a predictor of an adequate an- 
aesthetic dose. 

From international patent application no. WO 98/10701 (PCT/GB97/02435) 
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a control system and a method for calculating an index representation of 
the depth of anaesthesia are known. The method of calculating an index 
indication of anaesthetic depth is based on monitoring AEP produced by 
the patient and providing a signal corresponding to the coarseness of the 
5 monitored AEP signal. The raw AEP signal is divided into a series of 
sweeps or frames of a given duration, each sweep being synchronised with 
the repetitive audio stimuli. A number of sweeps n are recorded in se- 
quence and are averaged to produce a time average sweep. The anaes- 
thesia index is calculated for the time-averaged sweep. Each time a new 
10 series of sweeps is recorded, a new time-averaged sweep is determined 
from the most resent n sweeps, and the anaesthesia index for that time- 
averaged sweep is calculated. In this way the index is constantly updated. 

It has been observed that when a patient loses consciousness, the ampli- 
15 tudes of most AEP peaks are reduced and their latencies are generally also 
increased. These changes occur almost simultaneously and in the same 
direction with all patients. Therefore, a suitable index is one, which reflects 
these changes. 

20 An empirical algorithm has been developed for calculating the index, which 
algorithm is based upon the sum of the square roots of the difference be- 
tween every two successive points in the moving time-averaged sweep. 
This auditory evoked potential index is given by the following equation: 

255 

25 ^=*£VI*'-**t| 

w 

where Xi to x 2 56 are the sample points of the time-averaged frame and k is a 
scaling constant: 

30 The AEP index is calculated for every filtered time-averaged sweep, and a 
plot of the index against time can be generated for display and on a screen. 
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When the patient is awake, the index is typically in the range of 80 to 90, 
whereas during anaesthesia it is typically in the range of 35 to 40. 

An article in "Methods of Information in Medicine", 1996; 35: 256-260, with 
5 the title: "Autoregressive Modeling with Exogenous Input of Middle-Latency 
Auditory-Evoked Potentials to Measure Rapid Changes in Depth of Anes- 
thesia" by E. W. Jensen, P. Lindholm and S. W. Henneberg describe a sys- 
tem identification method, an autoregressive model with exogenous input 
(ARX), to produce a sweep-by-sweep estimate of the AEP. The method 

10 was clinically evaluated in 10 patients anaesthetized with alfentanil and 
propofol. The time interval between propofol induction and the time when 
the Na-Pa amplitude was decreased to 25% of the initial amplitude was 
measured. These measurements showed that ARX-estimated compared to 
MTA-estimated AEP was significantly faster in tracing transition from con- 

15 sciousness to unconsciousness during propofol induction (p < 0.05). 

It is the object of the invention to improve this measuring method in such a 
manner that a safer result is achieved considerably more rapidly, whereby 
the risk of treating a patient, e.g. by surgery, without full anaesthetization is 
20 reduced. 

The delay is reduced to about 6 seconds by using ARX modelling. 

It is a second object of the invention to make the procedure for anaestheti- 
25 zation more effective (time-efficient) and to reduce the staff workload. 

It is a third object of the invention to produce an apparatus for continuous 
monitoring of the level of consciousness, which apparatus is portable and 
easy to install and operate. 

30 

This is achieved by proceeding as disclosed in claim 1, e.g. using an appa- 
ratus as disclosed in claim 8. 
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Further advantageous characteristics are achieved by that which is dis- 
closed in the dependent claims. 

With the method and the apparatus according to the invention the possibil- 
5 ity is achieved of extracting the AEP with only a few repetitions, often as 
low as 15, which reduces the delay to approximately 6 seconds. As the 
AEP is a very complex signal comprising several peaks and troughs, it is 
desirable to map the AEP into a single number - an ARX index of easy in- 
terpretation, but containing the same information as the AEP. This is possi- 
10 ble by following the method and by using the apparatus according to the 
invention. The ARX index is typically larger than 60 when the patient is 
awake, and decreases when the patient is anaesthetized; a loss of con- 
sciousness will typically occur when the index gets below 28. 

15 The apparatus according to the invention can operate only with three sur- 
face electrodes, e.g. three propritary surface electrodes, the result being: 

- Fully updated AEP available within few seconds. 

- Significantly faster calculating AEP index than the traditional 
20 Moving Time Average method. 

- Significantly faster at tracing transition from unconscious- 
ness and vice versa. 

- Consistent, accurate readings. 

- Possible optimised display for use in operating theatre. 
25 - Touch screen - easy to operate and clean. 

- Fully graphic display. 

The method and the apparatus according to the invention monitor the level 
of consciousness during general anaesthesia independently of the biologi- 
30 cal variation of the patients with respect to tolerance and sensitivity of the 
anaesthetics. 



1 » 
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The method according to the invention can be performed on a computer as 
disclosed in claims 13 and 14. 

Brief description of the drawings 

5 

The invention will now be described in more detail with reference to the 
drawings, in which: 



Fig. 1 shows a schematic view of an embodiment of a system of the 
10 present invention for extracting the AEP (auditory evoked po- 

tentials) and for calculating an index, 

Fig. 2 shows a block diagram with the necessary procedures to pro- 
vide the extraction and index calculation of the AEP (auditory 
15 evoked potentials) according to an embodiment of the inven- 

tion, and 



Fig. 3 shows a block diagram illustrating the generalised ARX model 
(autoregressive with exogenous input), which is used accord- 
20 ing to the invention. 



Detailed description of the preferred embodiment 

In figure 1 is schematically shown a system 1 for extracting the AEP (audi- 
25 tory evoked potentials) and for calculating an index, which is indicative of 
the depth of anaesthesia of a patient. 

The patient 2 is subjected to repetitive sound signals, which are delivered to 
the patient by headphones, earphones or the like. These sound signals are 
30 in the form of "clicking" noise signals of short duration, approximately 1 - 2 
ms, which are delivered to both ears of the patient, and which produces 
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distinctive potentials, known as AEP or auditory evoked potentials in the 
electroencephalograph^ (EEG) response of the patient. The headphones 
and the equipment for producing the sound signals, e.g. a signal generator, 
are not shown in figure 1 . 

5 

Electrodes, with which EEG-signals (electroencephalographic signals) are 
provided, are attached to the head of the patient 2. Three or more elec- 
trodes, e.g. scalp electrodes, may be used. If three electrodes are used, 
they are attached to the patient 2 at the following positions: A positive elec- 
1 0 trode at the middle forehead, a reference electrode at the left forehead and 
a negative electrode behind the ear, preferably in the region of the mastoid 
process. Other positions may be used with similar results. 

The electrodes are connected to a patient cable, which leads to an amplifier 
15 3. This amplifier 3 is an instrumentation amplifier with a high common mode 
rejection ratio (CMMR). The amplified analogue signal is digitised by an 
A/D-converter (not shown) before being led to a digital processing unit 4. 
This digital unit may be in the form of a personal computer (PC), in particu- 
lar a single board personal computer. The digitised signal is analysed and 
20 stored on the computer 4, which is programmed according to algorithms for 
extracting the AEP (auditory evoked potentials) and according to algorithms 
for calculation of an index. These methods for extracting the AEP and for 
calculating the index shall later be described in further detail. 

25 The AEP and the index are displayed on a display 6, and with an input de- 
vice 5 is it possible to give instructions to the computer 4. This input device 
5 may be in the form of a touch plate or a touch screen, which can be com- 
bined with the display 6, in which way it will be possible to give instructions 
to the computer by touching the display. 

30 



The extracted and calculated values of the index and the AEP-signals can 
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be transmitted to external equipment (not shown) using a connection 7, 
which connection may be in the form of a RS-232 connection. The external 
equipment may for example be apparatus for administering drugs to the 
patient. Thus the system can be used for controlling the amount of drug 
5 being delivered to the patient in relation to the depth of anaesthesia. Other 
examples of external equipment are apparatus for monitoring or for treating 
the patient, for example apparatus for controlling the respiratory system. 

The procedures or functions included in a method and a system, called the 
10 A-line-monitor (AAI: A-line ARX Index), according to an embodiment of the 
invention is shown in the diagram in figure 2. 

The data representing the signals from the patient is introduced to the 
processing unit 4 in the form of a series of raw sweeps 11. In an embodi- 
15 ment of the invention these sweeps 11 are provided at a sampling fre- 
quency of 900 Hz, and each sweep consists of 70 samples, giving a sweep 
length of 80 ms. The click stimulus used for producing the signals was of 2 
ms duration and of an intensity of 70 dB above normal hearing level. 

20 Each sweep is initially processed by an artifact algorithm 12 in order to de- 
cide whether the sweep should be included for further processing. Two 
types of artifact algorithms are used. 

First, when artefacts are present in a sweep the amplitude of the raw data 
25 will in general be much larger than the amplitude of normal sweeps. The 
amplitude, expressed as a range of numbers, allowed of the A/D-converter, 
is for example 0 - 65534. The 95% EEG range is 15000 - 55000, hence if 
a sample is below 5000 or above 60000 a number of subsequent samples, 
e.g. 400 samples are rejected. 

30 

Second, during periods after saturation of the amplifier 3 the signal will 
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move in to a normal range, i,e. 15000 - 55000. However in this period a 
click artifact arising from the acoustic stimulus may be observed. Hence an 
algorithm that detects the click-artifact is implemented. The mean difference 
of the 10 samples before the click is calculated. If the difference between 
5 the first sample in a sweep n and the last sample in the sweep n-1 is larger 
than for example five times the previous calculated mean difference, hence 
a number of subsequent samples, e.g. the next 800 samples, are rejected. 

In order to improve the signal-to-noise ratio (SNR) a band pass filter 13 is 
10 included in the system previous to the application of the moving timer aver- 
age (MTA) and the autoregressive with exogenous input (ARX) model. The 
filter, which in an embodiment of the invention has a band pass from 16 to 
150 Hz, may be a fifth order Butterworth filter, which is digitally imple- 
mented. 

15 

The ARX model (autoregressive with exogenous input), with which a rapid 
extraction of the AEP is facilitated, shall now be described with reference to 
figure 3, which shows a generalised ARX model. 

20 The ARX model is obtained by adding an exogenous input to the AR model 
for analysing digital signals. Hence the ARX model is defined by the follow- 
ing equation: 

y(t) + ai • y(t- 1) + ... + a n ■ y(t- n) = bi ■ u(t) + + b m • u(t- m + 1) + e 

25 

where n is the order of the backwards coefficients (ai ... a n ) and m is the 
order of the forward coefficients (bi ... b m ). The output is y, u is the exoge- 
nous input and e is the error. 

30 On figure 3 is shown the AR-part 22, which is driven by white noise and is 
defined by the averaged EEG activity Pre-averaging may be done using 15 



WO 01/74248 



PCT/DK00/00636 



9 

sweeps. 

The exogenous input u to the block 23 is an AEP produced by averaging a 
number of sweeps, for example the latest 256 sweeps. The output y from 
5 the ARX model is an average of a number, preferably 15, of the latest col- 
lected sweeps consisting of averaged EEG background activity and AEP. 
When the coefficients of the model are determined, the ARX-AEP is ob- 
tained by IIR filtering of the exogenous input u. 

10 With a model order of for example five, the ARX equations are the following 

y6 = -aiy5 - .... -a 5 yi + biu 6 + + b 5 U2 + e 

y7 = -aiy 6 - .... -a 5 y2 + biu 7 + + b 5 u 3 + e 

15 

y 70 = -a^ - .... -a 5 y65 + biUyo + + b 5 U66 + e 

20 The error terms, e, are omitted and the equations are written on matrix 
form: 



25 



30 



-ys - y4 — y3 ->2 -y\ U6 us ua u3 u2 
-y6 -ys ->4 ->3 -yi ui us us ua U3 



->68 -J/67 -y66 -y6S U70 US9 U68 M67 Utej 



>7 



as 

bl 



WO 01/74248 PCT/DKOO/00636 

10 

The equation system shown above is an over determined set of linear 
equations. Hence Gaussean elimination or LU-decomposition will fail to 
give a satisfactory result. A very powerful way to solve an over determined 
system is the singular value decomposition (SVD), which solves the prob- 
5 lem in a least mean squares (LMS) sense. Singularities in a matrix, which 
means that the matrix does not have a full rank, often occur when the ma- 
trix is composed of data with no clear a priori knowledge. The singularity 
can occur if here is an ambiguity in the equations. This is a paradox be- 
cause on the one hand the system is over determined (more equations than 
10 unknowns) and at the same time it is undetermined because too many 
equations are linear (or close to linear) combinations of each other. SVD 
not only diagnoses and solves the singularity problems and produces a 
meaningful numerical result but it also provides the LMS solution. 

15 The model order is determined by considering the error function. The error 
function is defined by 

20 where e(i) = y(i) - y*(i) and N the number of samples in one sweep. The 
variables y(i) and y'Xi) are the real and the predicted pre-average of a num- 
ber of sweeps, for example 15 sweeps, respectively. The identification is 
tested by Andersons test on whiteness of the prediction error, e. If the pre- 
diction error is white at a confidence level of 95% the identification is ac- 

25 cepted. The optimal values. of n and m are selected by minimising the final 
prediction error (FPE) function defined by Akaike (Akaike H. "Statistical 
predictor identification", Ann. Inst. Statist. Math, 1970; No. 22, pp 203 - 
217): 



30 



FPE = L-(N + n + 1)/(N-n-1) 
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where n is the total number of coefficients of the ARX model. 

FPE represents a need of minimising the error function and the need of lim- 
iting the number of parameters of the ARX model. 

5 

The order of the ARX model should ideally be calculated for each sweep. 
This is a very time consuming process, hence to comply with the need of 
fast processing time an average model-order of five for both backward- and 
forward-coefficients was chosen. It is obvious that another suitable number 
1 0 may be chosen, preferably a number less than 1 0. 

As shown on figure 2, the output from the block 15, which is the AEP 
achieved by an moving time average (MTA) calculated over 256 sweeps 
with each sweep weighted equally, is led to the display 6 (fig. 1), where it 
15 may be displayed as an AEPivrrA-display 16. Further the output from the 
band pass filter 13, i.e. the EEG-signal, is led to the display 6 (fig. 1), where 
it may be displayed as an EEG-display 14. 

The output from the block 15 is also led to the ARX model 18 together with 
20 the output from the output from the block 17, which produces a MTA over a 
minor number of sweeps, preferably 15 sweeps. The output from the ARX 
model 18, which as explained above is an AEP extracted by the ARX and 
which is called an AEParx. is also led to the display unit 6 (fig. 1), where it 
may be displayed as an AEPARx-display 19. 

25 

In order to quantify the level of anaesthesia it is desirable to map the 
AEParx into an index. In order to do this, the output from the ARX model is 
led to a block 20 for index calculation, the function of which shall be ex- 
plained in the following. 

30 



The index calculated by the block 20 is called the A-line ARX-index (AAI) 
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and may be displayed on the display unit 6 (fig. 1), where it may be dis- 
played as an AAI-(AEP A Rx-)index display 21 . The AEP consists of several 
peaks. It is generally accepted that the amplitudes of the peaks with latency 
10-100 ms, corresponding to the middle latency AEP (MLAEP), de- 
5 creases when the patient is anaesthetized and at the same time the laten- 
cies of the peaks are prolonged. The index according to the invention pre- 
serves these two rules in order not to lose information. Furthermore in order 
to achieve a reliable index the following premises are complied with: 

10 1. Validity for the largest number of patients possible, independent of 

surgery and aesthetic drugs. 

2. Good dynamics is required between awake and asleep state in order 
to distinguish awareness changes from noise. 

15 

The AEP-index is calculated in a window of the AEP, which may preclude 
the start and the end of the window. Preferably the window is the 20 - 80 
ms window of the AEP and latency and amplitude changes in the AEP is 
weighted equally. The 20 ms start of the window is chosen not to include 

20 BAEP (Brainstem Auditory Evoked Potentials) and auricular muscular arte- 
facts and the 80 ms end of the window is chosen in order not to include 
LLAEP. This is because BAEP and LLAEP (Long Latency Auditory Evoked 
Potentials) do not correlate well to anaesthetic depth. The AEP index (AAl) 
according to the invention has shown good discrimination between con- 

25 scious and anaesthetized patients in previous studies. The AEP-index (AAl) 
is defined to reflect the hypnotic level during anaesthesia. 

First y is defined as: 
30 y=*i5>,-*/ + ,| 
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where X| are the samples of a sweep, 

ki are a constant, which preferably is larger than 0,0100 and less than 
0,02000, and which in particular may be chosen to be 0,0165, and 
where k2 and k3 are the start and the end samples, respectively, of the 
5 summation, chosen not to include the start and the end of the AEP window. 



If the AEP window consists of 70 samples, k 2 may preferably be 17, and k3 
may be 69. 



10 



The AAl-index is defined as: 



15 



20 



AAI -index = 



ify>60=>k A y+k s 
ify<l5andy>4=>k 6 +k 1 xe Y/p 



where k4, ks, ke, k 7 and p are constants. 

Preferably 0,2500 < k* < 0,3000, and in a most preferred form Id = 0,2786. 
Preferably 43,0000 < k 5 < 43,5000, and in a most preferred form k 5 : 
43,2857. 

Preferably 9,1000 < k6 < 9,8000, and in a most preferred form ke = 9,3769. 
Preferably 0,25 < k 7 < 0,30, and in a most preferred form k 7 = 0,28. 
Preferably 4 < p < 6, and in a most preferred form p = 5. 



25 In the most preferred embodiment of the invention the ARX index is defined 
as: 



0.0165 f>,-x j+1 | 



M7 



30 



AAI-index= j if ^>60=>0.2786>>+43.2857 

] ify<\5and y >4=*9.3169+0.2Sxe r ' ! 



M * 
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The index is in the range of 0 to 99, where increasing value indicates ele- 
vated level of consciousness. 

The method according to the invention may essentially be performed on a 
5 computer 4, e.g. a Single Board computer with 486 MHz clock frequency 
and provided with programming with matching software, e.g. using the pro- 
gramming language Borland Pascal, so that the method and the algorithms 
or parts thereof, which are explained in the description and disclosed in the 
method claims, are performed in such a manner that the desired result is ■ 
10 achieved, i.e. calculation and display of the AAl-index explained above, the 
program according to the invention is storable on any known type of com- 
puter readable medium, so that it is easily installed in the computer 4. 
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CLAIMS 



10 



15 



20 



25 



1. A method of extracting signals which are indicative of the level of con- 
sciousness of a patient, which method comprises the subjecting of the pa- 
tient to a repetitive audio stimulus, monitoring auditory evoked potentials 
(AEP) produced by the patient, extracting an AEP within a few repetitions, 
preferably more than 10 and less than 50, of the audio stimulus using an 
autoregressive model with exogenous input (ARX), and then calculating an 
index (AAI) indicative of anaesthetic depth. 

2. A method according to claim 1, characterized in that the method com- 
prises filtering of the AEP signals so that only signals with a frequency 
higher than 16 Hz and lower than 150 Hz are used in calculating the index. 

3. A method according to claims 1 and 2, characterized in that only signals 
in an AEP window consisting of about 100 ms (milliseconds) are used in 
calculating the index, preferably a window of 20-80 ms precluding the start 
and the end of the sweep. 

4. A method according to any of the claims 1-3, characterized in that the 
index is defined as: 



where k4, k 5 , k 6 , k 7 and p are constants, and where y is defined as: 
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where xi are the samples of a sweep, 

ki are a constant, which preferably is larger than 0,0100 and less than 
0,02000, and which in particular may be chosen to be 0,0165, and 
where k 2 and k3 are the start and the end samples, respectively, of the 
5 summation, chosen not to include the start and the end of the AEP window. 

5. A method according to claim 4, characterized in that the AEP window 
consists of maximum 80 samples, and in that k 2 may preferably be 17, and 
in that k 3 may preferably be 69. 

10 

6. A method according to claim 4 or 5, characterized in that k4, k 5 , k 6 , k 7 
and p are defined as follows: 

Preferably 0,2500 < k4 < 0,3000, and in a most preferred form k4 = 0,2786. 
15 Preferably 43,0000 < k 5 < 43,5000, and in a most preferred form k 5 = 
43,2857. 

Preferably 9,1000 < ke < 9,8000, and in a most preferred form ke = 9,3769. 
Preferably 0,25 < k 7 < 0,30, and in a most preferred form k 7 = 0,28. 
Preferably 4 < p < 6, and in a most preferred form p = 5. 

20 

7. A method according to any of claims 1-3, characterized in that the index 
is defined as: 



25 



AM - index = 



0.0165 Y,\x t -x M \ 



l /^>60=>0.2786y+43.2857 
ify<\5andy>4=>9.3769+0.28xe Y ' $ 



30 



8. An apparatus for extracting signals which are indicative of the level of 
consciousness of a patient comprising means for monitoring auditory 
evoked potentials (AEP) produced by the patient as a response to a repeti- 
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tive acoustic click stimulus, means for extracting an AEP within a few repe- 
titions, preferably more than 10 and less than 50 of the audio stimulus, 
means for using an autoregressive model with exogenous input (ARX), and 
means for calculating an index (Ml) indicative of anaesthetic depth. 

5 

9. An apparatus according to claim 8, characterized in that it further com- 
prises means for filtering the AEP signals, preferably a band pass filter with 
a band pass from 16 to 150 Hz, and preferably a digitally implemented But- 
terworth filter of fifth order. 

10 

10. An apparatus according to claim 8 or 9, characterized in that is com- 
prises means for providing sweeps at a sampling frequency of about 700- 
1000 Hz, preferably about 900 Hz, and an AEP window of maximum 100 
ms, in that only signals within the window are used in calculating the index, 

15 preferably a window of 20-80 ms precluding the start and the end of the 
sweep. 

11. An apparatus according to any of claims 8-10, characterized in that it 
comprises means for calculating an index (AAI) indicative of anaesthetic 

20 depth according to the methods described in claims 4-7. 

12. An apparatus according to claim 1 1 , characterized in that it comprises 
means to display or to signal the calculated index (AAI). 

25 13. A computer program comprising computer programming code means 
adapted to perform all the steps of any of the claims 1-7 when said program 
is run on a computer. 

14. A computer program as claimed in claim 13 and embodied on a com- 
30 puter readable medium. 
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